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■ , Abstract: Dark energy cosmologies with an equation of state parameter w less than —1 are 

Q I often found to violate the null energy condition and show unstable behaviour. A solution to 

' this problem may require the existence of a consistent effective theory that violates the null 

energy condition only momentarily and does not develop any instabilities or other patholog- 
ist \ ical features for a late time cosmology. A model which incorporates a dynamical scalar field 
\^ • ip coupled to the quadratic Riemann invariant of the Gauss-Bonnet form is a viable proposal. 



Such an effective theory is shown to admit nonsingular cosmological evolutions for a wide 
range of scalar-Gauss-Bonnet coupling. We discuss the conditions for which our model yields 
observationally supported spectra of scalar and tensor fluctuations, under cosmological per- 
turbations. The model can provide a reasonable explanation for the transition from matter 
dominance to dark energy regime and the late time cosmic acceleration, offering an interesting 
testing ground for investigations of the cosmological modified gravity. 
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1. Introduction and Overview 

Dark energy and its associated cosmic acceleration problem has been presented as three well 
posed conundrums in modern cosmology: (1) why the dark energy equation of state Wjy^ is 
so close to —1, (2) why the dark energy density is comparable to the matter density now, 
the so-called cosmic coincidence problem, and (3) why the cosmological vacuum energy is 
small and positive. Additionally, we do not yet have a good theoretical understanding of the 
result that the expansion of our universe is accelerating after a long period of deceleration as 
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inferred from the observations of type la supernovae, gravitational weak lensing and cosmic 
microwave background (CMB) anisotropies [1,2]; for review, see, e.g. [3]. 

A possible source of the late time acceleration is given by vacuum energy with a constant 
equation of state, known as the cosm,ological constant, or a variant of it, quintessence [4]. 
Several plausible cosmological scenarios have also been put forward to explain these problems 
within the context of modified theories of scalar-tensor gravity such as, coupled quintessence [5], 
k-inflation or dilatonic-ghost model [6], scalar-phantom model [7], Gauss-Bonnet dark en- 
ergy [8] and its various generalizations [9-11]. These ideas are interesting as, like quintessence, 
they offer a possible solution to the cosmic coincidence problem. The proposals in [7-11] are 
promising because they may lead to the observationally supported equation of state, w 
while more interestingly provide a natural link between cosmic acceleration and fundamental 
particle theories, such as superstring theory. In this paper we will focus primarily on the first 
two questions raised above within the context of a generalized theory of scalar-tensor gravity 
that includes the Gauss-Bonnet curvature invariant coupled to a scalar field if. 

The question of whether the gravitational vacuum energy is something other than a 
pure cosmological term will not be central to our discussion. But we note that Einstein's 
general relativity supplemented with a cosmological constant term does not appear to have 
any advantages over scalar-tensor gravity containing a standard scalar potential. The recent 
observation that the dark energy equation of state parameter tt; is ~ — 1 does not necessarily 
imply that the dark energy is in the form of a cosmological constant; it is quite plausible 
that after inflation the scalar field (p has almost frozen, so that T^(^) ~ const = A. In 
a cosmological background, there is no deep reason for expecting the energy density of the 
gravitational vacuum to be a constant, instead perhaps it can be determined by the underlying 
theory, as in the case where a scalar potential possesses many minima. It is thus worth 
exploring dynamical dark energy models, supporting both w < —1 and u) > — 1, and also 
w*{z) = dw/dz (where z is the redshift parameter). 

The expansion of the universe is perhaps best described by a monotonically decreasing 
Hubble expansion rate, implying that H = a/a — a? /a^ < 0, where a = a{t) is the scale 
factor of a four-dimensional Friedmann-Robertson- Walker universe and the dot represents a 
derivative with respect to the cosmic time t. In the presence of a barotropic fluid of pressure 
p and energy density p, this last inequality implies that the cosmic expansion obeys the null 
energy condition (NEC), p + p > 0, and hence w = p/p > —1. The standard view is that the 
Hubble expansion rate increases as we consider earlier epochs until it is of a similar order to 
that of the Planck mass, mpi ~ 10-*^^ GeV. However, in strong gravitational fields, such as, 
during inflation, the Einstein description of gravity is thought to break down and quantum 
gravity effects are expected to become important. This provides a basis for the assumption 
that the expansion of the universe is inseparable from the issue of the ultraviolet completion 
of gravity. In recent years, several proposals have been made in order to establish such a link. 
For instance, reference [12] introduced a system of a derivatively coupled scalar Lagrangian 
which violates the condition H < spontaneously: the model would involve a short-scale 
(quantum) instability associated with a super- luminal cosmic expansion (see also [13]). 

The beauty of Einstein's theory is in its simplicity. It has been remarkably successful as a 
classical theory of gravitational interactions from scales of millimeters through to kiloparsecs. 
Thus any modification of Einstein's theory, both at small and large distance scales, must be 
consistent with known tests. Several proposals in the literature, e.g. in [12,14], do not seem 
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to fall into this category as these ideas would involve modifications of Einstein's theory in a 
rather non-standard (and non-trivial) way. 

We motivate our work through the theoretical insights of superstring or M-theory as it 
would appear worthwhile to explore the cosmological implications of such models. In this pa- 
per we examine whether we can achieve observationally supported cosmological perturbations 
in the low-energy string effective action, which includes a non-trivial interaction between a 
dynamical scalar field ip and a Riemann invariant of the Gauss-Bonnet form, and study its 
phenomenological viability as a dark energy model. It is appreciated that a generalized theory 
of scalar-tensor gravity, featuring one or several scalar fields coupled to a spacetime curva- 
ture, or a Riemann curvature invariant, can easily account for an accelerated universe with 
quintessence {w > —1), cosmological constant {w = —1) or phantom [w < —1) equation of 
state without introducing the wrong sign on the scalar kinetic term. 

This paper is organised as follows. In the following section we discuss a general scalar- 
field Lagrangian framework and write equations of motion that describe gravity and a scalar 
field f, allowing non-trivial matter-scalar couplings. We discuss some astrophysical and cos- 
mological constraints on the model. In section 3 we present the construction of a number 
of scalar potentials for some specific constraints, in an attempt to gain insight into the be- 
haviour of the scalar potential for late time cosmology. In section 4 we discuss inflationary 
cosmologies for specific cases and study the parameters related to cosmological perturbations 
of the background solution. The problem of suitable initial conditions, given stable obser- 
vationally viable solutions with a full array of background fields for the general system is 
considered in section 5 using numerical and analytic techniques for both minimal and non- 
minimal scalar-matter couplings. In section 6 we present several remarks about the existence 
of superluminal propagation and/or short-scale negative instabilities for the tensor modes. 
Section 7 is devoted to the discussions of our main results. 



2. Essential Ingredients 

An unambiguous and natural way of modifying Einstein's General Relativity in four dimen- 
sions is to introduce one or more fundamental scalar fields and their interactions with the 

leading order curvature terms in the a' expansion, as arising in string or M theory compact- 
ifications from ten- or eleven-dimensions to four-dimensions. The simplest version of such 
scalar-tensor gravities, which is sufficiently general for explaining the cosmological evolution 
of our universe, may be given by 

S = Sgrav + Sjn, (2-1) 

with the following general action 

Viv) - ^/(^)^Gb) , (2.2) 

^{A\cp)ipm + Pr)), (2.3) 

where TZq^ = i?^ - 4 R^iy R'^'' + R^,ypa i?^'''*^ is the Gauss-Bonnet (GB) curvature invariant. 
We shall assume that 99 is a canonical field, so 7 > 0. The coupling fif) between ip and the 
GB term is a universal feature of all four-dimensional manifestations of heterotic superstring 



Sm = S{ip,A'^{(p)g^^,1prn) = 
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and M theory [15, 16]. For example, such a form arises at heterotic string tree-level if (p 
represents a dilaton, and at one-loop level if ip represents the average volume modulus; in a 
known example of heterotic string theory, one has /(v^) oc X]n=i e^""^)'^ > in the former 
case, while oc 99 — ^ e*^ + 4^^^^^ In (l — g-SnTre"^^ + ln2 < in the latter case (see also 
discussions in [17]). As discussed in [10, 11], a non-trivial or non-constant J{ip) is useful 
not only for modelling a late time cosmology but also desirable for embedding the model in a 
fundamental theory, such as superstring theory. Here we also note that from a model building 
point of view, only two of the functions V{ip), f{<p) and A((p) are independent; these can be 
related through the Klein-Gordan equation (see below). 

One can supplement the above action with other higher derivative terms, such as those 
proportional to /(</?) (Va'^V'^^?)^ and curvature terms, but in such cases it would only be 
possible to get approximate (asymptotic) solutions, so we limit ourselves to the above action. 
In the model previously studied by Antoniadis et al. [15], V{if) and Sm were set to zero. 
However, the states of string or M theory are known to include extended objects of various 
dimensionalities, known as 'branes, beside trapped fluxes and non-trivial cycles or geometries 
in the internal (Calabi-Yau) spaces. It is also natural to expect a non-vanishing potential 
to arise in the four-dimensional string theory action due to some non-pertubative effects of 
branes and fluxes. With supersymmetry broken, such a potential can have isolated minima 
with massive scalars, this then avoids the problem with runway behaviour of 99 after inflation. 

In this paper we fully extend the earlier analysis of the model in [10] by introducing 
background matter and radiation fields. We allow if to couple with both an ordinary dust- 
like matter and a relativistic fluid. We show that the model under consideration is sufficient 
to solve all major cosmological conundrums of concordance cosmology, notably the transition 
from matter dominance to a dark energy regime and the late time cosmic acceleration problem 
attributed to dark energy, satisfying ~ — 1. 

2.1 Basic equations 

In order to analyse the model we take a four-dimensional spacetime metric in standard 
Priedmann-Robertson- Walker (FRW) form: ds'^ = —dt'^ + a'^{t)Y^^=iidx^)'^, where a{t) is 
the scale factor of the universe. The equations of motion that describe gravity, the scalar 
field ip, and the background fluid (matter and radiation) are given by 

- ^ + 3^Hf^ ^ + + y + = 0, (2.4) 

4 (2e + 3) + - y - / - 2^i7/,^(l + e) + 'i>w^^^ = 0, (2.5) 
-f{p + 3Hp) + V^ + 3{l + e)l^H^-vQA''pm = 0, (2.6) 

where f = f^i^ipip'^ + f^,p(p, Vtp = dV/d(p, rj is a numerical parameter which we define below 
and 

H <P V{p) _pb PbA^ _ dlnAjp) 

where h stands for the background matter and radiation. Let us also define 

% ^ 3ff2 Q u = f^^H'^, Q.GB = 'pHf\^ = ux = p., (2.8) 
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so that the constraint equation ( p.4| ) reads Q^p + ^Igb + = 1- The density fraction Qh niay 
be spht into radiation and matter components: = 0,r + and Wb^b = Wm^m + Wr^r- 
Similarly, in the component form, Qpm = QiPm ■ Stiff matter for which Wm = 1 may also 
be included. The analysis of Steinhardt and Turok [18] neglects such a contribution, where 
only the ordinary pressureless dust {wm = 0) and radiation {wr = 1/3) were considered, in a 
model with f{if) = 0. The implicit assumption above is that matter couples to A'^{ip)g^i, with 
scale factor a, where a = aA{ip), rather than the Einstein metric g^,y alone, and pr oc 1/a^, 
so pr does not enter the (p equation of motion, (p.6|). That is, by construction, the coupling 
of (p to radiation is vanishing. This is, in fact, consistent with the fact that the quantity Q 
couples to the trace of the matter stress tensor, ^|^^r^^\ which is vanishing for the radiation 

component, = —p,^ + 3P^, = 0. We also note that, in general, p oc l/d^^^^'^\ thus, for 
ordinary matter or dust {w = 0), we have rj = 1, while for a highly relativistic matter {w = 1), 
we have ry = —2. Especially, for radiation (w = 1/3), we have 77 = 0; while, for any other 
relativistic matter having a pressure p > p > p/3, —2<r]<0. 

Equations (|2.4| )- (|2.6D may be supplemented with the equation of motion for a barotropic 
perfect fluid, which is given by 

^dpb I dpb p dpb 

where pb is the pressure of the fluid component with energy density pb- In the case of minimal 
scalar-matter coupling, the quantity dpb/dp ^ ^ and hence 

9.'f, + 2enb + ?,{l + Wb)^b = ^. (2.10) 

where the prime denotes a derivative with respect to AA^ = ln[o(t)] + const and a{t) is 
the scale factor of a four-dimensional Friedman-Robertson- Walker universe. In this case the 
dynamics in a homogeneous and isotropic FRW spacetime may be determined by specifying 
the field potential V{(p) and/or the scalar-GB coupling f{ip). 

For A{ip) 7^ const, there exists a new parameter, Q; more precisely. 

The variation in the energy densities of the ordinary matter, p^ , and the relativistic fluid, 

is) 

Pm , and their scalar couplings Qi are not essentially the same. Thus, hereafter, we denote 
the Q by Qd for an ordinary matter (or dust) and Qs for a relativistic matter (or stiff fluid 
coupled to if). Equation ( ^ ) may be written as 

Vl'^^^ + 2eVL^^^ + 3(1 + Wrr:)^m,s = -"illQ ,^m,s^' , (2.12) 

n'^^d + 2en^,d + 3nm,d = -^Qd^m,dV', (2.13) 
n'^ + 2enr + AQr = 0. (2.14) 



^Recently, Sami and Tsujikawa analysed the model numerically by considering Q ^ [19]. The authors 
also studied the case Q — const, but without modifying the Friedmann constraint equation. 
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For a slowly rolling scalar (or dark energy) field satisfying A99 oc AA^ (where = ln[a(t)] 
is the e- folding time), Qd is dynamically determined, for any /(</?), by the relation 

= QdV'^m,d- (2-15) 

For ~ ^, the scalar field presumably couples strongly to matter. However, for the 
present value of O^, ~ 10~^ (on large cosmological scales), one would hardly find a noticeable 
difference between the cases of non-minimal and minimal coupling of (p with matter. 
In the following we adopt the convention = SvrG = 1, unless shown explicitly. 



2.2 Cosmological and astrophysical constraints 

It is generally valid that scalar-tensor gravity models of the type considered in this paper 
are constrained by various cosmological and astrophysical observations, including the big 
bang nucleosynthesis bound on the (/^-component of the total energy density and the local 
gravity experiments. However, the constraints obtained in a standard cosmological setup 
(for instance, by analyzing the CMB data), which assumes general relativity, that is, (5 = 0, 
cannot be straightforwardly applied to the present model. The distance of the last scattering 
surface can be (slightly) modified if the universe at the stage of BBN contained an appreciable 
amount in energy in (^-component, like O,^ < 1/6 (see, e.g. [21]). 

The couplings Qd and Qs can be constrained by taking into account cosmological and 
solar system experiments. We should also note that observations are made of objects that 
can be classified as visible matter or dust [22,23]; hence a value of <C 1 is suggestive to 
agree with the current observational limits on deviations from the equivalence principle. If 
we also require <C 1, then the non- minimal coupling of Lp with a relativistic fluid may be 
completely neglected at present, since Jl^ <C ^d- In this subsection, as we are more interested 
on astrophysical and cosmological constraints of the model, we take fig = 0, without loss of 
generality. 

Under PPN approximations [24], the local gravity constraints on Qd and its derivative 
loosely imply that 

mliQl < 10-^ mpi \dQd/d^\ < 0{1). (2.16) 

If A{ip) is sufficiently fiat near the current value oi ip = (po, then these couplings have modest 
effects on large cosmological scales. Especially, in the case that A{(p) oc e^'^'^, the above 
constraints may be satisfied only for a small C (^ 1) since Qd = C is constant. On the other 
hand, if A{(p) oc cosh[(^Kp], then Qd is only approximately constant, at late times. Another 
possibility is that Qd oc e"''"'^; in this case, however, for the consistency of the model, one 
also requires a steep potential. In the particular case of A{p) oc e^'^^ with C ~ 0(1) (as one 
may expect from string theory or particle physics beyond the standard model), it may be 
difficult to satisfy the local gravity constraints (under PPN approximation), unless there is a 
mechanism like the one discussed in [25] (see also [26]). 

Only in the case Qdp' ^ 1; is the last term in (|2.13 ) possibly relevant, leading to a 



measurable effect including a 'fifth force causing violations of the equivalence principle (see, 
e.g., [22]). However, in our model, it is important to realize that the coupling Qd is dynami- 
cally determined by the ratio Qr/ip'^m)- Even if we allow a small slope for the time variation 
of the field, (p' ~ 0.01, which may be required for satisfying the BBN bound imposed on the 
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time variation of Newton's constant (see below), we find Qd ~ in the present universe, 
since ilr ~ 10"'^ and J7m ~ 0.3. In our numerical studies below we shall allow Qd in a large 
range 10~^ < < 0.01 and study the effect of Qd on the duration of matter (and radiation) 
dominance (cf figure ^). 

One may be interested to give an order of magnitude to the scalar field mass, rriip. This 
is especially because the lighter scalar fields appear to be more tightly constrained by local 
gravity experiments as they presumably couple strongly to matter fields. However, one of the 
important features of our model is that the dynamics of the field ip is not governed by the 
field potential V{ip) alone but rather by the effective potential (cf equation (^]B)): 



Here oc 1/a^, which is conserved with respect to the Einstein frame metric g^^- One can 
take (5 < 0, so that the last term above contributes positively to the effective potential; for 
the parametrization A{if) oc e'''^^"^^ , this requires C < 0. An important implication of the 
above potential is that the mass of the field ip, which is given by (= ^yVc{^, ipip), it not 
constant, rather it depends on the ambient matter density, pm- It may well be that, on large 
cosmological scales, ip is sufficiently light, i.e. rrifj, > 10~^^eV « (lO^^cm)"^, and its energy 
density evolves slowly over cosmological time-scales, such that Ves ~ V{4>) + | /('Z')^^ ~ 
10~^'^(GeV)^. However, in a gravitationally bound system, e.g. on Earth, where the pm is 
some 10'^'^ times larger than on cosmological scales, that is pm ^ Pc ^ 10^'^ ^ 2 g/cm^, the 
Compton wavelength of the field may be sufficiently small, Ac ~ m~ 1 < lO^eV"^ ~ 0.2 mm 
as to satisfy all major tests of gravity (see, e.g. [25], for similar discussions). 

There may also exist a constraint on the time variation of Newton's constant. With the 
approximation that the matter is well described by a pressureless (non-relativistic) perfect 
fluid with density pm, satisfying p^ <C Pm, which implies that <C 1 or A ~ const (cf 
equation 2.15| ), the growth of matter fluctuation may be expressed in a standard form: 



5 + 2(5i/ = 47rGefr 5, (2.17) 
where 5 = 6pm/pm- The effective Newton's constant Gefr is given by [27] 



Gcs — G 



l+3i}GB-^(t + ¥]^GB 



(2.18) 



where VLcb = Hf < 1. In the particular case of f{ip) oc e°"^, this yields 

Geff = G[1 + f{2Hip - (p)] . (2.19) 



It is not improbable that G^s — G for the present value of the field ipQ and the coupling 
fifo), since ip/ip <^ 1 and Hf < 1. Indeed, Ggg is suppressed (as compared to Ggs) by a 
factor of ip and thus one may satisfy the current bound '^^"^ / G^s < 0.01 //qj where Hq is 
the present value of H, by taking ip/H < 0.01. This translates to the constraint [Gnow — 
Gnucieo I /Gnow (^now " ^nucieo) < 10~^^yr~^. Another opportunity for our model to overcome 
constraints coming from GR is to have a coupling f{ipo) which is nearly at its minimum. 
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The constraint on the time variation of Newton's constant may arise even if i^cB = 0, 
given that Qd is nonzero. In this case we find 



Ges = GA\^){l + Ql). 
Nevertheless, in the present universe, since pr <C Pm, ^iv^) ~ const and is smah (<C 1). 

3. Construction of inflationary potentials 

In this section we study the model in the absence of background matter (and radiation). 
Equations ( |2.4| )-( |2^ ) form a system of two independent equations of motions, as given by 

7x2 + 2y + 6/i-6 = 0, (3.1) 
/i' + (e - 1)^ - 7a;2 - 2e = 0. (3.2) 

The solution e = corresponds to a cosmological constant term, for which w = —1 . We 
shall assume that the scalar potential is non-negative, so y > 0. Evidently, with > 0, as is 
the case for a canonical (p, the inequality p, < 1 holds at all times. Here /U = 1 is a saddle point 
for any value of e. Thus, an apparent presence of ghost states (or short-scale instabilities) at 
a semi-classical level, as discussed in [28] , which further extends the discussion in [29] , is not 
physical. We shall return to this point in more detail, in our latter discussions on small scale 
instabilities incorporating superluminal propagations or ghost states. 

An interesting question to ask is: can one use the modified Gauss-Bonnet theory for 
explaining inflation in the distant past? In search of an answer to this novel question, it is 
sufficiently illustrative to focus first on one simple example, by making the ansatz 

x = xoe"^/2^ p = poe'^^, (3.3) 

where xq, pq, a and /3 are all arbitrary constants. The implicit assumption is that both x 
and p decrease exponentially with A^, or the expansion of the universe. We then find 

From this we can see that a transition between the e > and e < phases is possible if 
p 0- Moreover, the first assumption in ( |3.3| ) implies that aN = 2lnip + lnipi, where ipi is 
an arbitrary constant. The scalar potential is then given by 

Vi^) = Hi^f (3 - Xoip' - Ai (^2/3/^) , (3.5) 

where Ao,Ai are (arbitrary) constants. For a canonical ip, so > 0, we have Aq > 0, 
whereas the sign of Ai is determined by the sign of pQ. Note that a = /3 is a special case 
for which the potential takes the form V{ip) = H'^i^p) (3 — (Aq + Xi)(p'^). The quantity e 



■^The equation of state (EoS) parameter 111 = ^ — — ^^,^^'2^ — — 1 — § £• The universe accelerates for 



w < —1/3, or equivalently, for e > —1. In fact, for the action (2.2), it is possible to attain w < —1 or 



equivalently e > 0, without introducing a wrong sign to the kinetic term. 
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(and hence 1 + w) cannot change its sign in this case. Typically, if /3 = 2a, then the scalar 
potential would involve a term which is fourth power in if, i.e., V oc i/^((/?) cp^. At this point 
we also note that the potential ( |3.5D is different from a symmetry breaking type potential 



F oc (A lb m^(/9^ + Xiif'^ + • • •) considered in the literature; here it is multiplied by H'^{ip). An 
inflationary potential of the form 1/ oc c/?^ is already ruled out by recent WMAP results, at 
3fj-level, see e.g. [34]. In the view of this result, rather than the monomial potentials, namely 
V{if) oc {ip/ipoY , a scalar potential of the form V{ip) oc H'^{ip){(p/ipoy, with p > 2, as implied 
by the symmetries of Einstein's equations, is worth studying in the context of inflationary 
paradigm. Further details will appear elsewhere. 

Another interesting question to ask is: is it possible to use the modified Gauss-Bonnet 
theory to explain the ongoing accelerated expansion of the universe? Before answering this 
question, we note that, especially, at late times, the rolling of p can be modest. In turn, it is 
reasonable to suppose that ip/H const, or x ~ xq. Hence 

V{p) = ^ (e - W + 6/io e'^*) = H{pf (Ao + Ai e^-^) , (3.6) 
where (j) = {ip — 'Pq)/xq and the Hubble parameter is given by 

/ ^ ,N (2+a;o)/2/3-l 



with Hq being an integration constant. Interestingly, a non- vanishing /((/?) not only supports 
a quartic term in the potential, propotional to H'^{p) 93^, but its presence in the effective 
action also allows the possibility that the equation of state parameter w switches its value 
between the w > —1 and w < —1 phases. We shall analyse the model with the choice ( |3.6| ) 
and in the presence of matter fields, where we will observe that the universe can smoothly 
pass from a stage of matter dominance to dark energy dominance. 

In the case where (/? is rolling with a constant velocity, ip' = c, satisfying the power-law 
expansion a{t) oc t^/^, or equivalently H = HQe~^^ and p ^ 1, we find that the scalar 
potential and the scalar-GB coupling evolve as 

V = Voe- 2P<^ + 3(3p + l)fiH^ e (3.8) 

f = he^^'^-^e^'+'P^* + f2, (3.9) 
-"0 

where cp = {p — po)/c, fi and /2 are arbitrary constants, and 

_ {6-6p + 5c^-pc^)Hi , _ 2p-c^ 

2(1 +p) 2p{p + l)H^ 

The potential is a sum of two exponential terms; such a potential may arise, for example, 
from a time-dependent compactification of 10 or lid supergravities on factor spaces [30]. In 
general, both V{p) and f{p) pick up additional terms in the presence of matter fields, but 
they may retain similar structures. In fact, various special or critical solutions discussed in 
the literature, for example [8], correspond to the choice /i = /2 = 0. 

We can be more specific here. Let us consider the following ansatz [8]: 



aoc(t + ti)", p = P0 + pin{t + ti), (3.11) 



- 9 - 



for which obviously both e and are constants: 



e^H/H^ = -l/a, ^' ^ ^ = ^. (3.12) 

H a 

For e < 0, one can take ti = 0; the Hubble parameter is given by 

H = \a\e'^^'P^'~^-^''\ (3.13) 
The scalar potential is double exponential, which is given by 

y = 1/oe-2<^ + yie("-^)'^, (3.14) 

where (p = (ip — 9?o)//3 and 

Vo = TT/ — T-TT ) H = -3(a + 3)aci. (3.15) 

2[a + 1) 

The scalar-GB coupling /((/?) may be given by 

/ = /oe2^ + /ie(-+3)^ + /2, (3.16) 

with 

^°=2(l + a)a2' ^^ = a2(a + 3)' = ^^"^^^ 

Of course, the numerical coefficient /2 does not contribute to Einstein's equations in four 
dimensions but, if non-zero, it will generate a non-trivial term for the potential. 
In the case e > (or H > 0), the above ansatz may be modified as 

aoc(too-i)", ip = ipo + (3ln{too-t), (3.18) 

where a < 0. The Hubble parameter is then given by -fT = — a(too — i)^^- Such a solution 
to dark energy is problematic. To see the problem, first note that although this solution 
avoids the initial singularity at i = 0, it nonetheless develops a big-rip type singularity at an 
asymptotic future t = t^o- This is not a physically appealing case. Also the above critical 
solution may be unstable under (inhomogeneous) cosmological perturbations, which often 
leads to a super-luminal expansion and also violates all energy conditions. 

The reconstruction scheme presented, for example, in ref. [31] was partly based on some 
special ansatz, e.g. ( ^.11]) , which may therefore suffer from one or more future singularities. 
However, as we show below, for the model under consideration there exists a more general 
class of cosmological solutions without any (cosmological) singularity. 



4. Inflationary cosmology: scalar-field dynamics 

Inflation is now a well established paradigm of a consistent cosmology, which is strongly 
supported by recent WMAP data ^. It is also generally believed that the small density (or 

^The small non-uniformities observed for primordial density or temperature fluctuations in the cosmic 
microwave background provide support for the concept of inflation. 
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scalar) fluctuations developed during inflation naturally lead to the formation of galaxies, stars 
and all the other structure in the present universe. It is therefore interesting to consider the 
possibility of achieving observationally supported cosmological perturbations in low-energy 
string effective actions. For the model under consideration, this can be done by using the 
standard method of studying the tensor, vector, and scalar modes. We omit the details of 
our calculations because they are essentially contained in [29] (see also [28,32]). 

4.1 Inflationary parameters 

One may define the slow roll parameters, such as and ijh, associated with cosmological 
perturbations in a FRW background, using two apparently different versions of the slow roll 
expansion. The first (and more widely used) scheme in the literature takes the potential 
as the fundamental quantity, while the second scheme takes the Hubble parameter as the 
fundamental quantity. A real advantage of the second approach [33] is that it also applies to 
models where inflation results from the term(s) other than the scalar field potential. Let us 
define these variables in terms of the Hubble parameters ^: 

(in the units k = 1). Here, as before, the prime denotes a derivative w.r. to e-folding time 
N = ln[a(t)], not the field ip. One also defines the parameter ^h, which is second order in 
slow roll expansion: 

«„4(^^^)"^(«,„-V2^^)"\ (4.2) 

These definitions are applicable to both cases, V{ip) ^ and F((/?) = 0, and are based on the 
fact that inflation occurs as long as ^(^) < holds. The above quantities (e/f , 77^ j Cif) are 
known as, respectively, the slope, curvature and jerk parameters. 

Typically, in the case /(</?) = 0, or simply when ir^Gsl so that the coupled GB 

term becomes only subdominant to the field potential, the spectral indices of scalar and tensor 
perturbations to the second order in slow roll expansion may be given by [33] (see also [34]) 

nn-\ = -4eH + 2riH - 2(1 + c)e% - ^(3 - 5c)eHVH + ^(3 - c)^|, 

riT = -2eH-iS + c)el + {l + c)eHVH, (4.3) 

where c = 4(ln 2 + 7) — 5 0.08. For the solutions satisfying — 0, implying that both 
eH and tjh are much smaller than unity (at least, near the end of inflation) and their time 
derivatives are negligible, we have 

n-ji-'i-- -4:eH + 2riH, nr ^ -^en- (4.4) 

In fact, n-ji and ny, along with the scalar-tensor ratio r, which is given by r « 16eH + 32c(eff — 
flH)^Hi are the quantities directly linked to inflationary cosmology. 



*The slow roll variable is tn, not e, the latter is defined by e = H' /H. 
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Below we write down the results directly applicable to the theories of scalar-tensor gravity 
with the Gauss-Bonnet term. Let us introduce the following quantities: ^ 



where 



H (p x' + ex q' 



Even in the absence of the GB coupling = 0), hence £4 = 65 = eg = 0, there are particular 
difficulties in evaluating the spectral indices n-fi and ut, and the tensor-scalar ratio r, in full 
generality. In perturbation theory it is possible to get analytic results only by making one or 
more simplifying assumptions. In the simplest case, one treats the parameters almost as 
constants, so their time derivatives are (negligibly) small as compared to other terms in the 
slow roll expansion. An ideal situation like this is possible if inflation occurred entirely due to 
the power-law expansion, a(r) oc |r|~ "^/^^'''^^ where the conformal time r = —l/[aH{l + e)]. 
In this case, the spectral indices of scalar and tensor perturbations are well approximated by 



n-Te - 1 = 3 



3 + ei + 262 + 2e4 
1-ei 



riT 



3 - ei + 2e6 
1-ei 



(4.7) 



Note that not all are smaller than unity. Nevertheless, an interesting observation is that, 
for various explicit solutions found in this paper, the quantity €4 is close to zero, while €5, 
eg can have small variations during the early phase of inflation. But, after a few e-folds of 
inflation, AA^ > 5, these all become much smaller than unity, so only the first two terms (ei. 



€2) enter into any expressions of interest. In any case, below we will apply the formulae (4.7) 
to some explicit cosmological solutions. 

The above relations are trustworthy only in the limit where the speeds of propagation 
for scalar and tensor modes, which may be given by 

2 _ 1 , ^^[4e(l-/x) + /-^] 2 _ 1-/ _ 1 - + g/^ 8^ 

+(l-M)[27(l-/.)x2 + 3^2]' ^^-^7/^- l-n ' ^^-^^ 

where ^ = fH and x = (f/H, take approximately constant values. These formulae may be 
expressed in terms of the quantity z^, which is defined below in ( 4.1C| ), using the relation 



^ = v' — 2ev. The propagation speeds depend on the scalar potential only implicitly, i.e. 
through the background solutions which may be different for the V{'-p) = and V{^) 7^ 
cases. In the case where and are varying considerably, the derivative terms like e's, €5 
are non-negligible, for which there would be non-trivial corrections to the above formulae for 
riTi and nx- Furthermore, the spectral indices diverge for ei ~ 1, thus the results would apply 
only to the case where ei <C 1. In this rather special case, which may hold after a few e-folds 



^The parameter 63 defined in [29] (see also [32]) is zero in our case since the action (2.1) is already written 
in the Einstein frame, and does not couple with the Ricci-scalar term in this frame. 
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of power-law expansion through to near the end of inflation, we find that the tensor-to-scalar 
ratio is approximately given by 

..|:«16?24Lz4±V/£eV. (4.9) 



This is also the quantity directly linked to observations, other than the spectral indices n-jz 
and riT- The WMAP data put the constraint r < 0.55 at 2a level. 

4.2 Inflating without the potential, V{ip) = 

In this subsection we show that in our model it is possible to obtain an inflationary solution 
even in the absence of the scalar potential. It is convenient to introduce a new variable: 

i^iip) ^ f{v)H\ (4.10) 

This quantity should normally decrease with the expansion of the universe, so that all higher- 
order corrections to Einstein's theory become only sub-leading ^. In the absence of matter 
fields, and with V{ip) = 0, the equations of motion reduce to the form 

^, ^_6(l + e) + (5 + e)7x2 



27X 

u' = l + 2eu--^x'^. (4.11) 
6 

Fixing the (functional) form of i^(v) would alone give a desired evolution for e. Notice that for 
X = const, one has e = const; such a solution cannot support a transition between acceleration 
and deceleration phases. It is therefore reasonable to suppose that x ^ const, and also that 
1/ varies with the number of e-folds. A'". Take 

oc e'"^, (4.12) 

where m is arbitrary. As the simplest case, suppose that m ~ or equivalently u ~ const 
= i/Q. In this case we find an explicit exact solution which is given by 

e = -/3o + Ptanh{(3AN), (4.13) 

where AA' = A -|- Aq, with Nq being an integration constant '', and 

fi^irJ^^-S-^-'o + l (4.14) 

Clearly, the case f = must be treated separately. The evolution of (p is given by ip' = 
±-^(6 -I- 12efo)/7. The Hubble parameter is given by 

H = Ho e"*^ cosh(/3AA). (4.15) 



^Particularly in the case u ~ const, the coupled Gauss-Bonnet term \f{^)'Tix3B = •if{v)H'^{l + e) = 
3i>oH^{i^ + e) may be varying being proportional to the Einstein-Hilbert term, R/2 = 3H^{2 + e). 

'^One may choose the constant No differently. Here we assume that the scale factor before inflation is oo, 
so initially N = ln(a(t)/ao) = 0. As the universe expands, iV S> 0, since a{t) ^ ao- This last assumption 
may be reversed, especially, when one wants to use the model for studying a late time cosmology, where one 
normalizes o(t) such that its present value is ao, which implies that iV < in the past. 
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For uq > 0, both Po and /3 are positive and Po > P, implying that Hubble rate would decrease 
with the number of e-folds. In turn, the EOS parameter w is greater than —1. However, if 
i/Q < 0, then we get /3 < 0, whereas Po < for -0.2 < uq < and Po > for uq < -0.2. The 
vq < solution clearly supports a super-luminal expansion, see also [35]. 

With the solution ( p^ , we can easily obtain 50 — 60 e-folds of inflation (as is required to 
solve the horizon and the flatness problem) by choosing parameters such that 1 + /3 — /3q > 0. 
An apparent drawback of this solution is that it lacks a mechanism for ending inflation. We 
will suggest below a scenario where this problem may be overcome. 




Figure 1: The index n-ji — 1 as the function of i^o, with V{(p) = (cf equation 4.13); uti does not vary 
much with the number of e-folds except that it is blue-tilted for AA^ < 0. Even though f o < generally 
implies a super-luminal inflation, it is possible to achieve a red-tilted spectrum for i>q > —2/3. 



The constraint on the spectral indices might give a stringent bound on the scalar-Gauss- 
Bonnet coupling f{ip), or the GB density fraction, i^cB- The scalar index satisfies the bounds 
provided by the recent WMAP results, namely n-ji = 0.958^g'Q}g, for both signs of but in 
a limited range, see Fig. In the V{ip) = case, in order to end infiation the function 
z^((^) could be modified in a suitable way; in the case mN < 0, implying that z^(<^) decreases 
with N, infiation ends naturally. With this particular example, we do not wish to imply that 
the model with a vanishing potential alone is compatible with recent data. Such a model 
may suffer from one or more (quantum) instabilities at short distance scales and/or local GR 
constraints [27] once the model is coupled to matter fields. As in [36], we find that a model 
with V{(p) > is phenomenologically more viable than with V{ip) = 0. 

As noted above, a super-luminal expansion may occur if < 0, which corresponds to the 
case 5 < studied before in [15]. However, unlike claimed there, we find that it is possible to 
get a cosmic inflation of an arbitrary magnitude without violating the (null energy) condition 
H/H^ < 0. For the solution (4.13), this simply required vq > 0; the Hubble rate would 



decrease with N. For uq < 0, one gets H/H^ > 0, especially, when AN > 0. The precise 
value of redshift z (= e^^ — 1) with such a drastic change in the background evolution actually 
depends on the strength of the coupling f{^), or the time variation of ip. 
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4.3 Inflating with an exponential potential 

Let us take the scalar potential of the form 

= Vbe-2c^/'^o, (4.16) 
leaving /(v) unspecified. The system of autonomous equations is then described by 

dN ifQ 



dx _ {I + e){y - 3) 2y (5 + e)x 
dN 7x 7(/3o 2 ' 

1 y 7X 

u = . 

X 3x 6 

These equations admit the following de Sitter (fixed point) solution 



(4.18) 
(4.19) 



x = 0, V = Ao, H = ^^, u = u{y^), (4.20) 

which corresponds to the case of a cosmological constant term, for which w = —1. In the 
particular case that n oc x, so that u = f^ipH"^ is varying (decreasing) being proportional to 
the kinetic term for if, we find an explicit solution, which is given by 

. = ^o-a^otanha(iV + Ar^), n = -f, , = ^l^'^o -2xi6 + ,x^) ^ 

3 (po{6 + jx'^) 



where a = — < 1 ^^'^ is an integration constant. This solution supports a 

transition from deceleration (e < —1) to acceleration (e > —1), but inflation has no exit. The 
model should be supplemented by an additional mechanism, allowing the field (/? to exit from 
infiation after 50 — 60 e- folds of expansion (see a discussion below). 

For a negligibly small kinetic term, x ~ 0, we have ilcB oc ~ 0, while for a slowly 
rolling If, i^GB 7^ 0. This behaviour is seen also from our numerical solutions in the next 
section. The evolution of (p is given by 

ip = J xdN = Nipo - (fio lncosha(iV + N^) + const. (4.22) 

The scale factor, as given by a cx exp[J' H dt] oc exp [J" exp[J edN] dt\ , is regular everywhere. 
The coupling /((/?) = f[u{^) /ip'H"^] dip is rather a complicated function of ip, which is generally 
a sum of exponential terms. But given the fact that x ~ const at late times, or when 
\a{N + A''^)! > 2, we get /^^ oc e^'^^'^°. It is not difficult to see that the special solution 
considered in [8] corresponds to this particular case. To see it concretely, take 

/((^) = /oe2'^/'^o, (4.23) 



or, equivalently ipou' — 2u{ipoe + x) = 0. For the solution ( 4.21| ) this last condition imposes 
the constraint = 6 and hence 

X = xq, u = — 7x0/3. (4.24) 
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The choice made in [8] for V{(p) and /((/?) thus over-constrains the system unnecessarily. 

Can the solution ( 4.21D be used for explaining early inflation? In order to answer this 
question concretely, one may wish to evaluate cosmological parameters such as n-ji and tit- 
We can get the value of n-ji in the range 0.94 < n-ji < 0.99 by taking 9 < ^/^'^o ^ 20. In 
fact, riTi may exceed unity precisely when e > 0, leading to a momentary violation of the 
null energy condition (cf figure ^ (left panel)). However, after a certain number of e- folds of 
expansion, the tensor modes or the gravitational waves damp away (cf figure § (right panel)), 
leading to the observationally consistent result |nr| < and rin — 1^0- leading to slightly 
red-tilted spectra (cf figure 0). 
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Figure 2: The left panel shows that the equation of state parameter w that momentarily crosses 
the value w — —1. From top to bottom — 4,5,6,8 and 10. The right panel shows the tensor- 

to-scalar ratio of cosmological perturbations; to adhere to the experimental constraint coming from 
recent WMAP results, namely r < 0.5, one might require an almost flat potential, y/y(fQ > 16. 




Figure 3: The spectral indices n-ji — 1 and riT evaluated for the solution (4.21) by using the formulae 
given in (4.7); from top to bottom ^(fo = 10,8,6 and 4. 

For the validity of ( [4. 7] ) , we require not only e'j <C 1 but also that the speed of propagation 
for scalars and tensor modes, c| and c^, are nearly constant. This is however not the case 
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for the above solution when < AA^ < 3, for which both and change rapidly. That 
means, the relations ( ^.7] ) must be modified in the regime where c^,c^ are varying. After a 
few number of e-folds of expansion, we have c^, ~ const. In this limit the spectral indices 
are found to be slightly red-tilted, i.e. n-ji — 1 < 0, for ^/^'^PQ > 8. 

A few remarks may be relevant before we proceed. The WMAP constraints on inflationary 
parameters that we referred to above, such as, the scalar index n-ji and the tensor-to-scalar 
ratio r, are rather model dependent and may differ, for instance, for a model where the 
"inflaton" field is non-minimally coupled to gravity, or the spacetime curvature tensor. Hence, 
although those values may be taken as useful references, they may not be strictly applicable 
to all inflationary solutions driven by one or more scalar fields. In fact, some of our cases may 
have more freedom (less tight constraints) than the canonical models with a scalar potential 
proportional to or rn^ip^, and /(v^) = 0. We leave for future work the implications of 
WMAP data to scalar potentials of the form V oc H^{(p)(p" (n > 2). 

4.4 Inflating with an exponential coupling 

Let us take the scalar-GB coupling of the following form: 

/,^ = /oe2^/^«, (4.25) 

but without specifying the potential, V{(p). With this choice, the system of monotonous 
equations is given by 

^ 2(x + e^ (4.26) 

|| = ?i±I£! + (l_3.).-^. (4.27) 

dl\ u (fQ 

y = 3--x^ - 3ux. (4.28) 



The de Sitter fixed point solution for which x = is the same as (|420| ). In the following we 
consider two special cases: 

(1) Suppose that fi = const, that is, 0,gb = const = fiQ. Then we find 

^ /ip + dx ^ 2{-fipox^ - (2 - fJ'o)x^ + 2i^oipox) 

^ ~ //o - 2 ' dAT ~ (2 - fio)ipo ' ^ ' ' 

These equations may be solved analytically for //q = 2/3 ^; as the solutions are still messy 
to write, we only show the behavior of w in Mathematica plots. In the next section we will 
numerically solve the field equations, in the presence of matter, allowing us to consider all 
values of /xq- We should at least mention that the above system of equations has a pair of 
fixed point solutions: 

2 - Mo - \/(2 - Mo)^ - 2 - /^o + \/(2 - /^o)^ - 87V?o/Wo „n\ 
xi = , X2 = . (4.30) 

2-jipo 2jipo 



*This is actually a critical point in the phase space, which may be seen also in cosmological perturbation 
analysis, see, e.g., [20]. 
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The fixed point xi is an attractor, while the X2 is a repeller provided that xi,X2 are real. In 
the case where Sjipo^J-o > (2 — /io)^ the solution diverges. The solution also diverges for initial 
values of x such that x < and x > X2- If these conditions are not violated, the solutions 
would always converge to the attractor fixed point x = xi. The evolution of the solution is 
monotonic from the initial value of e to the e given by the attractor fixed point at x = xi. 
The initial value of w for a wide range of /io and initial values of x may be read from figure 
§(a), and evolve to the w given in figure §(b) for a specific value of /io- 
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Figure 4: The left panel (a) shows the initial equation of state parameter w with various initial 
conditions for flcB = fJ-o f^nd x, while, the right panel (b) shows a possible variation of w for a 
particular value of /^o! solutions evolve to these values regardless of Xini provided the fixed points exist 
and < Xini < X2- We have taken ip'^^ = 2 and 7 = 1 

(2) Next suppose that u = const = uq, instead of ux = fi =const. In this case, the 
quantity x decreases quickly with the expansion of the universe; the GB density fraction, 
^GB, also decreases with the number of e-folds. The explicit solution is 

e = -^ x = (431) 

(/3o' 7^30 + 1*0 + (2 - uo¥'o)iii e(2/«0'P0-l)Af ' 

where ui is an integration constant. Without any surprise, as x — > 0, e — > 0; the model then 
corresponds to the cosmological constant case. 

Inflation is apparently future eternal for most solutions discussed above, e.g. ( 4.13| ) and 



( ^.21| ). However, in a more viable cosmological scenario, the contribution of matter field 
may not be completely negligible. This is because, during infiation, the field 0, while slowly 
evolving down its potential, which satisfies Acp cx AA^, decays into lighter particles and 
radiation. The inflaton may even decay to heavier particles, especially, if the reheating of the 
universe was due to an instance preheating [37]. In turn, a significant amount of the energy 
density in the (p component might be transforming into the radiation and (baryonic plus dark) 
matter, with several hundreds of degrees of freedom and with all components present, e.g. 
stiff matter (p = p) and radiation (3p = p)- In turn, the slow roll type parameter e would 
receive a non-trivial contribution from the matter fields. Explicitly, we find 

e = -^{l + w)n,-^, (4.32) 
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where ^Ib = + ^^^^ ^ = (^'m ~^Pr)/^Pm + /^r)- Inflation ends when there is significant 
portion of matter fields, which makes e more and more negative and hence e < — 1. 

Before proceeding to next section we also wish to make a clear separation between the 
inflationary solutions that we discussed above and the dark energy cosmologies that we will 
discuss in the following sections. Although it might be interesting to provide a natural link 
between the early universe inflation and the late time cosmic acceleration attributed to dark 
energy, postulating a string-inspired model of quintessential inflation, the time and energy 
scales involved in the gravitational dynamics may vastly differ. As we have seen, through the 
construction of potentials, the scalar potential driving an inflationary phase at early epoch 
and a second weak inflationary episode at late times could be due to a single exponential term 
or a sum of exponential terms, but the slopes of the potential of the leading terms could be 
very different. In fact, one of the very interesting features of an exponential potential is that 
the cosmological evolution puts stringent constraints, e.g. during big-bang nucleosynthesis, 
on the slope of the potential, but its coefficient may not be tightly constrained. That is, even 
if we use the potential V{ip) = Vq e~^^, with /3 ~ 0(1 — 10), for explaining both the early and 
the late time cosmic accelerations, the magnitude of the coefficient Vq can be significantly 
different, hence indicating completely different time and energy scales. 



5. When matter fields are present 

The above results were found in the absence of radiation and matter fields. It is thus natural 
to ask what happens in a more realistic situation, at late times, when both radiation and 
matter evolve together with the field ip. In such cases a number of new and interesting 
background evolutions might be possible. Also some of the pathological features, like the 
appearance of super-luminal scalar modes, may be absent due to non-trivial scalar- matter 
couplings. 

5.1 Minimally coupled scalar field 

We shall first consider the case of a minimal coupling between (p and the matter, so A{ip) = 1. 



Then, equations (|23)-(|p|) and (|2lc| ) can be easily expressed as a system of three independent 



equations. After a simple calculation, we find 

6 = -^(lnl^,)'-^(l + u;6), (5.1) 

2 

^' = 3(1 + Wkfib) + 2e + ^- y-(2 + e)/i, (5.2) 

subject to the constraint + 0,gb + = 1- Note that the background equation of state 
parameter Wb is a function of A^, not a constant. One can write Wb^b = Wm^ln^ + WmQm + 
Wr^rj so that Wm = and Wr = 1/3 for pressureless dust and radiation, respectively, while 

(s) 

for a relativistic fluid, e.g. hot and warm dark matter, we have Wm > 0. In fact, the first 



equation above, (5J), gives the condition for accelerated expansion, which requires e > — 1. 
There are two ways to interpret this equation: either it gives in terms of Wb, assuming 
that the ratio Q^Qb is known, or it provides for flxed Wb the value of that may be related 
to the ratio QyClb- Of course, in the absence of kinetic term for ip [so x = 0), f{ip) is uniquely 
fixed once V{(p) is chosen. Such a construction is however not very physical. 



- 19 - 



An interesting question to ask is: What is the advantage of introducing a non-trivial 
couphng f{f), if it only modifies the potential or the variable y, without modifying the 
functional relation between e and 0,b? For the model under consideration without such a 
coupling there does not exist a cosmological solution crossing the phantom divide, w = —1. 

As a very good late time approximation, suppose that the ratio ip/H takes an almost 
constant value, say P, rather than 99 = 0, along with H cx e~'^'^^'^, where (p = \/l{f~fo)/P = 
A'^. The background matter density fraction is then given by 

^Ih = ^0 exp[A0 — 3(1 + Wb)4'], (5-3) 

where CIq = c\ {6wb + 8 — A)/[3(2 — A)]. For A < 2, the physical choice is ci > 0, so that 
Ofe > 0. The scalar potential and the GB density fraction evolve as 

^ ^ 3(2-A)+/3^(5-A/2) ^_ _ ^_ 3(1+^^)^ ^ ^ ^^^^^ 
2 -|- A 

2 + A 2 — A 3 

where ci,C2 are integration constants. The second term in each expression indicates a non- 
trivial contribution of the background matter and radiation. The presence of ci term in 
( |5.4D implies that the value of the density fraction $7^ must be known rather precisely for an 
accurate reconstruction of the potential. Equation ( ^.51) implies that 



^ /o + /i e2^<^-3(i+-i>)'^ - /2 e(2+V2)</.^ (5.6) 

where /i and /2 take, respectively, the signs of ci and 02- V is sufficiently flat near the 
current value of if, then both the background and Gauss-Bonnet density fractions, fi^ and 
^GB, may take nearly constant values. At late times, Wb — Wm — for the pressureless dust, 
while the radiation contribution may be neglected since VLr ~ 10~^. 

5.2 Non-minimally coupled scalar field 

To investigate a possible post-inflation scenario, where the scalar field may couple non- 
minimally to a relativistic fluid or stiff-matter other than to ordinary matter or dust, we 
consider the following three different epochs: (1) background domination by a stiff relativistic 
fluid, (2) radiation domination and (3) a relatively long period of dust-like matter domination 
which occurred just before the current epoch of dark energy or scalar field domination. To 
this end, we define the fractional densities as follows: Vlg for stiff-matter, for radiation and 
for dust, where the respective equation of state parameters are given hy Wg = 1, Wr = 1/3 
and = 0. Equations ( |2.12| )-( ]2?T^ ) may be written as 

fi; + 2e^s + 317s(l + Ws) = -3r],Qsnsip', (5.7) 
+ leQd + 30rf = -3Qdndip'. (5.8) 
0; + leflr + = 0, (5.9) 

In order for current experimental limits on verification of the equivalence principle to be 
satisfied, the coupling must be small, <C 1, at present, see e.g. [22]. 
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Superstring theory in its Hagedorn phase (a hot gas of strings), and also some brane 
models, naturally predict a universe filled with radiation and stiff-matter. To this end, it 
is not unreasonable to expect a non-trivial coupling of the field (p with the stiff-matter: 
highly relativistic fluids may have strong couplings of the order of unity, Qs ~ 0{1). More 
specifically, as we do not observe a highly relativistic stiff-matter at the present time, i.e. 

~ 0, the scalar-stiff-matter coupling Qg in the order of unity is not ruled out and may 
have a quantifiable effect on the evolution of the early universe. 

5.3 The ^GB = const solution 

It is interesting to study the case of a constant ^gb as it allows us to evolve the potential 
without constraining it through an ansatz. An exponential potential is generally used in the 
literature due to the simplification it affords allowing a change of variables and hence a system 
of autonomous equations. 

It is worthwhile to investigate whether a physically interesting result can be found without 
the need for an ansatz for the potential. To this end, let us take ^gb — const = fiQ. With 
the ansatz ( [4.25| ), we have 



potential V = yH^ and the effective potential A{ip) = V{ip) + 3f{ip)H'^{l + e) explicitly. We 
can see that there is a strong correlation between the length of dust-like matter domination 
and the parameter values, no and 0,d^im- 

The dust domination period is expected to last for about 5 — 7 e-folds of expansion. 
However, in the above case, no parameter ranges meet this condition while still maintaining a 
(relatively) long period of radiation domination. This is problematic as at least some period 
of radiation domination appears to be a requirement for reconciling any cosmological model 
with observations. Our results here are therefore presented in the spirit of a toy model, which 
allow the evolution of an unconstrained potential for parameter ranges and give solutions 
with the expected qualitative features. 

Figure ^ shows the period of dust domination observed for various values of ^o- We see 
that, in the parameter range /io ^ 10^'^, there exist solutions supporting both radiation and 
stiff-matter domination prior to matter domination. While for larger values of /Uq the Gauss- 
Bonnet contribution suppresses the background fields and we may observe a quick transition 
to scalar-field domination, missing either the radiation dominated era or the matter dominated 
era, or both. For smaller values of /Uq, no solutions exist for the parameters used above due 
to constraints in the field equations. The behaviour of Wcs for varying /xq and the general 
evolution of the cosmology are shown in figure ^ and figure 0, respectively. In the cases where 




(5.10) 





where 
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Figure 5: The period of dust-like matter domination as a function of fj,o for f^^ — f^e^^, 7 = 1, 
Q2 = iQ-b g^j^j g2 ^ 0.01. The initial conditions are x = 10"'^, y = 10-l^ Vld = 6.1 and fi^ = 0.45. 




Figure 6: The behaviour of Wcs for varying ^0 for /.^ = foe^'^, 7 = 1, = 10 ^ and Ql = 0.01. 
The initial conditions are x = 10^^, y = 10^^^, ild = 0.1 and ilr = 0.45. 

the second derivative of Wes is monotonic there is no radiation domination {wes ~ 1/3) or 
dust domination (wjefr ~ 0). For lower values of ^q, the second derivative of obviously 
changes sign and hence implies a period of background matter or radiation domination. We 
show the behaviour of the scalar potential and the corresponding effective potential in figure ^. 

5.4 Simplest exponential potentials 

We now wish to consider the evolution of the full system while putting minimal restrictions on 
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Figure 7: The evolution of the fractional densities and Wcff (solid line, black) for a fixed /io = 10~^ 
with parameters 'y = 1, ipo = 2/3, = 10^^ and Ql = 0.01, and initial conditions as a; = 10~^, 
y = 10^"'^^, fid = 0.1 (dots, brown) and fir = 0.45 (dashes, green), fig is represented by long dashes 
(orange) and il^p by dot-dash (red). 




Figure 8: The left panel shows the evolution of V{(p) while the right panel shows the evolution of 
Veff = V{ip) + 3f{(p)H^{l + e) for no = 10-^ 7=1, f{(p) cx e^-^, Ql = lO^^ and ^ q.OI with initial 
conditions x = IQ-'^ , y = 10"l^ Qd = 0.1, Qr = 0.45, H = 0.01 and fiip)^^^ = 10"^". 

the evolution of the cosmological constituents. To do this we employ simple single exponential 
terms for both the field potential and the scalar-GB coupling: 

/,^(¥.) =/oe"'^ and y((^) = Vb e"^'^. (5.13) 

The commonly invoked exponential potential ansatz has some physical motivation in super- 
gravity and superstring theories as it could arise due to some nonperturbative effects, such 
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as gaugino condensation and instantons. One can see that these choices are almost undoubt- 
edly too naive to allow all the expected physical features of our universe from inflation to 
the present day. This is because generally the slopes of the potential considered in post- 
inflation scenarios are too steep to allow the required number of e-folds of inflation in the 
early universe. As a post-inflation approximation, however, these may hold some validity, 
as one can replicate many observable physical features from nucleosynthesis to the present 
epoch while allowing non-trivial scalar-matter couplings. In follow-up work we will discuss 
the possibility of a two-scalar fields model where the potential related to one scalar field meets 
the requirements of inflation while the other scalar drives the late time cosmology. 

The ansatze (5.13) allow us to write an autonomous system 



dx 1 

TT7 = -TT [27a;e + 67X - 2/3y + 6^(1 + e) - 67?,Q,f], 
aiV z7 



dy_ 

dN 
du 

dN 



+VdQd{7x'^ + 21/ + 6{ux + nr + ns- 1)}] (5.14) 
-y(/3x + 2e) (5.15) 

u{ax + 2e) (5.16) 



where 



6{2'yux — 27 — 

+6"fux — 672/ + 37^x^ + 18u^ 
-VdQdC^l^x^ + 18u - ISu'^x - Guy 

- {rjdQd - risQs) (I27X - 18uO 



I87 (1 + Wd — Wd^r — Wd^s + Wr^r + Ws^s) 



o 2 2 
87 WdX 



Gf3uy — Gjaux"^ — ISjWdUx — GjuudV 
12jx + Ajxy + 12'jxftr + 9jux'^) 

(5.17) 



which along with the continuity equations ( |5.7D -( ^l9| ) and the Friedmann constraint equation, 
+ fi^ + i^GB + f^r + = 1) allow us to proceed with numerical computation. 

We look particularly at two values for /3 here, /3 = \/2/3 and /3 = \/3, while keeping the 
other parameters constant to limit the vast parameter space. These values may be motivated 
by various schemes of string or M theory compactifications (see, e.g. [30]). We will discuss 
the effects of varying these other parameters and some of the physically relevant results. The 
evolution of the various constituents is shown in figures |g|-|lO|. 

The /? = \/2/3 case appears to show a smoother evolution and has a short period during 
which there is acceleration and an appreciable amount of matter in the universe. In the table 
below we look at some of the features of the accelerating period for the solution given in 
figure |9[ At the onset of acceleration, Wes = —1/3, we have Qd = 0.650. Within the best-fit 
concordance cosmology the present observational data seem to require Weff < —0.74 [39]. 
However, in the above case, this requires less matter (0^ < 0.12) than in ACDM model. 
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Figure 9: Evolution of the fractional densities and Wcfi (solid line, black) where (} = •\/2/3, 7 = 1, 
a = 12, Ql = 10"^ and Ql = 0.01 with initial conditions a; = %/6 x IQ-*, y = 5 x IQ-^o, u = 0.08, 
fir — 0.549 (dashes, green) and Vlg = 0.45 (long dashes, orange), fid is represented by dots (brown), 
r2(p by dot-dash (red) and ^gb by dot-dot-dot-dash (blue). 



Condition imposed at present 


^-^accelerating 


Implied w^q 


Implied 


Qrf = 0.27 


0.582 


-0.648 


N/A 


Weff = -0.74 


0.955 


N/A 


0.123 


-Weff = - 0.9 


5.48 


N/A 


3.02 X 10"^ 


WeS = - 1.0 


5.55 


N/A 


2.46 X 10"^ 


AiV = 0.69 


N/A 


-0.685 


0.215 


AiV = 0.91 


N/A 


-0.734 


0.134 



The recent type la supernovae observations [40] appear to indicate that the universe may 
be accelerating out to a redshift of z ~ 0.4 — 1 [41]. In terms of the number of e-folds of 
expansion, this corresponds to IS.N = ln(l + z) ^ 0.34 — 0.69, if one assumes a rescaling of 
^ at the present epoch using the freedom in choosing the initial value of the scale factor, 
oq. We have chosen the initial value of y, or the ratio V{cp)/H'^, such that the period of 
dust-like matter domination is ^Nq^ ~ 6.5 as this corresponds to a total redshift of z ~ 1100 
to the epoch of matter-radiation equality. 

The evolution with f3 = \/3 does not seem to be in favourable agreement with current 
cosmological observations. The solution may undergo a period of sudden change in when 
the Gauss-Bonnet contribution becomes appreciable (or significant); this would have to occur 
around at the present epoch as to retrieve the current value of u^e// ~ — 1- This result does 
not reconcile with either the supernovae data, which seem to indicate a longer period of 
acceleration, or with a constraint for the present value of Gauss-Bonnet density [27] which 
may require that i^cB ^ 0.2, or with the presently anticipated value of (~ 0.27). 
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Figure 10: Evolution of the fractional densities and Weff (solid line, black) where (3 — \/3, 7 = 1, 
a = 12, Ql = lO"'"^ and ^ q.OI with initial conditions x = V6 x lO""*, y = 5 x lO^^o^ ^ ^ q.OS, 
ilr = 0.549 (dashes, green) and fig = 0.45 (long dashes, orange), ild is represented by dots (brown) 
and by dot-dash (red), floB by dot-dot-dot-dash (blue). 
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Figure 11: Evolution of WcS with /3 = 1/2/8, 7 = 1, ^ IQ-^ and = 0.01; a = 12 (dash-dot, 
brown), 8 (dashes, green), 3 (dots, red), y'2/3 (solid, black), with initial values x — \/6 x 10"^, 
y^5x 10-20^ u = 0.08, fir = 0.549 and fl, = 0.45. 

We can observe an oscillatory crossing of w^fi = — 1 limit for all cases in which the 
Gauss-Bonnet contribution becomes appreciable, even momentarily. Such a behaviour may 
be seen in variants of scalar-tensor models [42]. In our case, the amplitude of these oscillations 




Figure 12: Evolution of Woff witli /5 = ^3, 7 = 1, ^ 10"^ and = 0.01; a = 12 (dash-dot, 
brown), 8 (dashes, green), 3 (dots, red), (solid, black), with initial values x — a/G x 10~^, 

?/ = 5 X 10^20^ y ^ 0.08, ^ 0.549 and = 0.45. 




Figure 13: The period of dust-like matter domination for varying and Qd for 7 = 1, a = 12 and 
/3 = 1/2/3 with initial conditions x = ^6 x 10-^ y = 5 x IQ-^o, u = 0.08, = 0.549 and = 0.45. 



corresponds to the amplitude of the oscillations seen in the Gauss-Bonnet contribution and 
hence is heavily dependent on the slope of the scalar-GB coupling; for large a we observe 
much larger oscillations. These oscillations damp quickly as the Gauss-Bonnet contribution 
becomes negligibly small, and settle to a late time evolution for which Wgfj « — 1. As this 
limit is approached from above, none of the issues inherent with super-inflation or a violation 
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Figure 14: Evolution of Woff for different values of Q^: = 10^ (dots, green) 100 (dashes, orange), 
(solid, black) and 1 (dash-dot, red) with fixed Qd = VlO'^ and7 = l,a = 12,/3 = ^/2J3. The 
initial conditions are x = x 10^'', y = 5 x 10^^", u — 0.1, fir — 0.549 and Qs = 0.45. The larger, 
more violent oscillations seen for Qg > 1 are due to significant scalar-stiff-matter couplings, while the 
smooth oscillations for smaller couplings at some later stage are due to an appreciable contribution of 
the coupled Gauss-Bonnet term. 

of unitarity will be applicable to the late time cosmology. 

As the parameter space for the initial conditions is very large we have presented solutions 
with initial conditions and parameters selected to give reasonable periods of radiation and 
dust-like matter domination as well as other physically favourable features. Although the 
quantitative behaviour is observed to change smoothly with changes in initial conditions and 
parameters, such changes do have some qualitative effects as limits of certain behaviour are 
encountered. A quantitative variation in the period of dust-like matter domination can be 
attributed to altering the initial value of y or r^^^ini; lower values of yini extend the period 
before scalar field domination begins. This is effectively a change in the initial potential and 
has a monotonic effect on the epoch at which scalar field domination begins. In a qualitative 
sense, it has no effect on the period of radiation domination until a ymi is selected which 
is large enough that the scalar field contribution completely suppresses the dust-like matter 
domination period. The value of ^d,im is relevant to the epoch of matter-radiation equality, 
larger values result in an earlier epoch. For small values, no dust-like matter domination 
occurs and the solution is entirely dominated by the other four constituents considered. The 
limits at which all these effects occur are dependent on other parameters and hence discussion 
of actual values instead of the general behaviour does not add further insight. 

We note that the exponential terms for both the potential and the coupling have some 
noticeable effects on the evolution of the system. This may be seen to an extent in terms of 
P in figures It does, however, appear that the ratio of these parameters also influences 

the expansion during dust-like matter domination remaining, relatively unchanged when this 
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Figure 15: Variation of WcS for non-zero Qd as compared to the Qd = case, Awoff = Wcs{Qd 
0) — Wcff{Qd = 0); = 10~^ (solid, black), 10^^ (dots, green) and 10~* (dashes, orange) and 
Qs = 0.1. Clearly, for higher values of the scalar dust-matter couplings, > 10~^, we find slightly 
larger variations in Weff ■ 



ratio is constant within realistic /? and a parameter ranges. Phenomenological bounds on 
these values have been studied in [20]. 

For a ~ /3 we do not observe a significant period of Gauss-Bonnet contribution and hence 
no crossing of the Wgfr = ~1 hmit. The closer the value of a is to this limit the later this 
period of significant Gauss-Bonnet density fraction occurs. As a increases there is a minimum 
epoch at which the Gauss-Bonnet contribution becomes significant, this epoch occurs after the 
scalar field becomes the dominant component in the energy budget of our universe. Results 



for various values of a are shown in figures 11 and 12 



The consequences of matter coupling to the scalar field on the cosmology are of particular 
interest. The effect in the period of dust-like matter-domination is minimal in relation to 
both and Qs and can be seen in figure 13, though we allowed the value Qs ^ Qd- As 
there is a lack of observational or theoretical motivation refuting the possibility of a high 
relativistic-matter-scalar coupling we consider a range that extends beyond Qs = 1 whereas 
the dust-like-matter-scalar coupling must take values < 10~^ due to the current level of 
experimental verification in solar system tests of GR. 

The solution undergoes a transition in qualitative behaviour when we extend through 
the limit Qs,iim ~ where the exact onset and amplitude of its effects are dependent on 
the other parameter values taken. Hence our remarks again apply to the general behaviour 
observed rather than specific cases. Couplings greater than Qs.Um cause a non-negligible 
reemergence of the stiff-matter contribution at the end of dust-like matter domination. The 
density fraction of the relativistic matter undergoes a damped oscillations, generally with a 
much shorter period of oscillation than the oscillations observed in the non-negligible Gauss- 
Bonnet contribution. This causes corresponding oscillations in the effective EoS, while still 
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generally showing the same overall trend as for lower Qg values of the otherwise same solution, 
with the density fraction of the relativistic matter stabilising to a non-zero late time value. 
This cosmological behaviour does not appear to be physically valid as no mechanism to 
generate this relativistic matter seems plausible and would hence lead us to suggest that the 
value of Qs would have an upper-bound such that Qs,max ~ Qs,iim- The effects of Qs and Qd 
on the effective equation of state parameter, if eff, may be seen in figures 14 and [l^. Note 
that in figure 15 we have considered Aw^ff rather than w^ff as the effects of varying 
are so minimal that no discernable variation can be seen otherwise. In these figures we have 
considered the same solutions as in figure ^. The magnitude of the effects are, however, the 
same throughout the parameter space. 




Figure 16: Evolution of the fractional densities and WcS (solid line, black) for the potential ansatz 
V{if) = iJ2(Ao + Aie"'^'^) where 7 = 1, a = 9, /3 = ^/2/3, Aq 10-^ Ql = 10"^ and = 0.01 with 
initial conditions x = V6 x 10^^, y ~ 9x 10^^, u = 0.1, il^ = 0.549 and = 10^^. is represented 
by dashes (green) , by dots (brown) by long dashes (orange) , fl^ by dot-dash (red) and flc b by 
dot-dot-dot-dash (blue). 



5.5 A canonical potential 

Motivated by our results in section 3, we consider the potential 

V{^) = H\Ao + Aie'^'P). (5.18) 



Imposing this ansatz along with the scalar-GB coupling given in equation ( 5.13| ), /(</?) oc e' 



aO!<fi 



allows us to find solutions that have reasonable agreement with concurrent observations while 
not crossing the Weg = — 1 limit at any stage of the evolution. From figure (|^) we can see that 
the cosmic evolution shows a smooth progression to Wes = — which may be physically more 
sensible. The amplitude of the Gauss-Bonnet density fraction at maximum is dependent on 
the values of slope a, for smaller a it never becomes relevant. The period of dust-like matter 
domination again shows a heavy dependence on the initial value of y or the potential V{ip). 
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Here we also make an essential remark related to all our numerical plots in sections 5 
and 6. Taking into account the relation AA^ = ln(l + z), one may (incorrectly) find that the 
present epoch corresponds to AA^ ~ 15. However, in our numerical analysis, we are interested 
in the general evolution of the universe, in terms of the e-folding time A^ = ln[a(t)], from the 
(early) epoch of radiation-dominated universe to the present epoch. Indeed, our solutions are 
completely invariant under a constant shift in AA^, so the choice of Nq in AA^ = A^ + Nq is 
arbitrary. One may choose Nq such that ~ 0.27 at present, or alternatively ~ 15, so 
that = 1 (the present value of the scale factor) corresponds to A^ ~ 0. 

6. Remarks on ghost conditions 




Figure 17: The speeds of propagation for scalar and tensor modes, corresponding to the solution 
(4.21) with ^y^^po = 3,4,5 and 6 from top to bottom (bottom to top) for the left (right) panel. 

Recently, a number of authors [28, 38] ^ discussed constraints on the field-dependent 
Gauss-Bonnet couplings with a single scalar field, so as to avoid the short-scale instabilities 
or superluminal propagation of scalar and tensor modes These conditions are nothing but 

< < 1, < 4 < 1, 1 - ^ > 0. (6.1) 

As long as < 1 and < e, we also get < 1 and < 1. The short-scale instabilities 
observed in [28] corresponding to the value /i ~ 1, or equivalently u' — 2eu ~ 1, may not be 
physical as the condition /i ~ 1 effectively invalidates the assumptions of linear perturbation 
theory (see e.g. [20,29]). Furthermore, in the limit — > 1, other higher order curvature 
corrections, like cubic terms in the Riemann tensor, may be relevant. To put this another 

®The model studied in [38] may find only limited applications within our model, as the kinetic term for ^ 
was dropped there, and also a rather atypical coupling /(i^s) oc tp^ , with n < 0, was considered. 

^'^If the existence of a superluminal propagation is seen only as a transient effect, then the model is phe- 
nomenologically viable. However, a violation of causality or null energy condition should not be seen as an 
effect of higher curvature couplings, at least, for late time cosmology; such corrections to Einstein's theory are 
normally suppressed as compared to the Ricci-scalar term, as well as the scalar potential, at late times. For 
a realistic cosmology the quantity f{(p)H^ should be a monotonically decreasing function of the number of 
e-folds, iV, since f{ip)RGB oc f{ip)H^ and V{ip) oc H'^{ip). 
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Figure 18: The variation of witli AA^; (top panel) P = -y/2/3 and a — 

12 (dasli — dot, brown), 8 (dashes, green), 3 (dots, orange), \f2j% (solid, black); (bottom panel) /3 — \/3 
and a = 12 (dash — dot, brown), 8 (dashes, green), 3 (dots, orange), (solid, black). Other param- 

eters are chosen as 7 = 1, = 10"^ and = 0.01, with initial conditions x = -\/6 x 10"**, 
y = 5 X 10^20, u = 0.08, Vtr = 0.549 and 17,, = 0.45. 



way, the Gauss-Bonnet modification of Einstein's theory alone is possibly not sufficient for 
describing cosmology in the regime ^gb ^^totai- For the consistency of the model under 
cosmological perturbations, the condition = IOgbI < 2/3 must hold, in general [20,27]. 

Let us first neglect the contribution of matter fields and consider the special solution 
( [4.2l| ), obtained with a simple exponential potential. From figure (17) we see that both 
and temporarily exceed unity, because of which cosmological perturbations may exhibit a 
superluminal scalar or tensor mode. However, this behaviour can be significantly different in 
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Figure 19: The variation of with AiV; (top panel) f3 = ■y/2/3 and a = 

12 (dash — dot, brown), 8 (dashes, green), 3 (dots, orange), ^^2/3 (sohd, black); (bottom panel) /3 — \/3 
and a = 12 (dash — dot, brown), 8 (dashes, green), 3 (dots, orange), -^^2/3 (solid, black), for 7=1, 
Ql = 10-^ and Ql = 0.01 with initial conditions a; = ^6 x 10"'', y = 5x lO'^o, u = 0.08, fir = 0.549 
and Qs = 0.45. 

the presence of matter, especially, if ip is allowed to couple non-minimally to matter fields. 

The propagation speeds normally depend on the spatial (intrinsic) curvature of the uni- 
verse rather than on a specific realization of the background evolution during the stage of 
quantum generation of scalar and tensor modes (or gravity waves). Thus a plausible expla- 
nation for the existence of super-luminal scalar or tensor modes {c^ > 1 or > 1) is that 
at the initial phase of inflation the spatial curvature K is non-negligible, while the interpre- 
tation of c|i (and c^) as the propagation speed is valid only for K = 0. Additionally, in the 
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presence of matter, the above result for does not quite hold since the scalar modes are 
naturally coupled to the matter sector. The result for may be applicable as tensor modes 
are generally not coupled to matter fields. 

Prom figures (jl^) and (19) we can see that, for ansatze such as ( ^.13 ), and may 
become negative, though temporarily, if one allows a larger slope for the scalar-GB coupling, 
namely a ^> /3. It is precisely this last case for which and may also take values larger 
than unity at subsequent stages, leading to super luminal propagation speeds for scalar and/or 
tensor modes (see also [19,43] for discussions on a similar theme). The case where > 1 
normally corresponds to the epoch where the contribution of the coupled Gauss-Bonnet term 
becomes significant. However, for smaller values of a and /3, satisfying a < f5, both and 
never become negative and there do not arise any ghost-like states, though depending upon 
the initial conditions the tensors modes may become superluminal, temporarily. 



7. Discussion and Conclusions 



Considering earlier publications on the subject by one of us (IPN), perhaps it would be useful 
first to outline the new results in this paper. Originally [10] the model was analysed by 
making an appropriate ansatz for the effective potential, namely 

V{if) + ^ f{^)RlB = A((^) ^ H\^) [3 + e((^)] , (7.1) 

which is a kind of gauge choice, and is well motivated by at least two physical assumptions: 
firstly, the effective equation of state < 1 for e > —3 secondly, this choice would allow 
us to find an explicit solution and compute the primordial spectra more easily. In this work 
we have relaxed such a constraint and extended our analysis of the cosmological solutions 
of these systems by allowing non-trivial couplings between the scalar field and the matter 
fields. The cosmological viability of such a generalized theory of scalar-tensor gravity is fully 
investigated by placing minimal constraints on the model parameters and the scalar-matter 
couplings. 

In summary, we outline the main findings of the paper. We discussed some astrophysical 
and cosmological constraints applicable to a general scalar-tensor gravity models and then 
outlined the conditions under which those constraints may be satisfied for the model under 
study, with generic choices of the model parameters. Under the assumptions that the quantity 
</?' = (p/H and the time-derivative of the scalar-GB coupling decrease during inflation being 
proportional to some exponential functions of = ln[a(t)], namely <p/H oc e"i^ and fH oc 
e"2^, the (reconstructed) scalar potential was shown to take an interesting form V{(p) = 
H'^{lp){Cq + C-^ip^ + C2 9j"2/"i). Such a potential, being proportional to the square of the 
Hubble parameter, would naturally relax its value as the universe expands and might have 
useful implications for the early universe inflation. With the approximation y?' ~ const, V{'p>) 
was shown to take the form V{'~p) = H'^{(p) (Vq + e^'^), for which the effective equation of 
state approaches —1 without exhibiting any pathological features at late times. 

We have studied the cases where inflaton field (p rolls with a constant velocity (99' = c) 
and/or the universe experiences a power law inflation, a{t) oc a°^. In several cases of physical 

^^Only the stiff-matter can saturate this hmit, so A{ip) = 0, or equivalently e{ip) — —3. 
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interest, the form of the (reconstructed) potentials and the scalar-Gauss-Bonnet couplings 
are found to be a sum of exponential terms, rather than a single exponential. We have 
explicitly shown that, for a spatially flat metric K = 0, non-singular inflationary solutions 
are possible for a wide range of scalar-GB couplings. We have found explicit cosmological 
solutions by imposing minimal restrictions on model parameters; for example, we allowed the 
scalar potential to take the simplest exponential form, but without specifying the coupling 
The model was also shown to admit an epoch of inflation for a vanishing potential, 
V{ip) = 0. This should not come as a surprise since we have already known from Starobinsky's 
work [44] that inflation can easily be triggered by the 4D anomaly driven i?^-terms. However, 
in this case, although the spectral indices of cosmological perturbations may be brought close 
to observationally supported values, such as 0.96, it is generally difficult to obtain an 

inflationary solution supporting a small tensor-to-scalar ratio, r <^ 0.5. 

We have shown (both analytically and numerically) that getting only a small deviation 
from w = —1 in our model necessarily implies a non-trivial scalar-GB coupling. The question 
of course is how natural is the existence of the phase w < —1 and the phase supporting matter 
to dark energy dominance. Depending upon the initial values of the model parameters, a 
momentary crossing of the EoS w = —1 may be realized in various contexts. We have 
explained under what conditions u; -|- 1 can be negative, given the required period of matter- 
dominance is maintained before realizing the currently observed acceleration of the universe. 
We have also established rather systematically that in theories of scalar-tensor gravity with 
the coupled Gauss-Bonnet term, the quantity w + 1 can be negative only momentarily, while 
approaching the limit w = —1 from above. Hence indicating that none of the issues inherent 
with a super-luminal propagation or a violation of unitarity will be applicable for a late time 
cosmology. It would be of interest to extend our analysis to the K = ±1 cases. 

Our analysis is not conclusive in answering whether the model under consideration is 
viable for describing our universe from the epoch of nucleosynthesis until today without any 
violation of unitarity, for instance, due to an imaginary propagation speed for the tensor or 
scalar modes. Such violations imply the existence of some sort of ghost states, or superluminal 
propagation, implying a break down of the weak energy condition or causality. Nevertheless, 
such pathological features of the model observed for simplest choices of the potential, e.g. 

(X e°"^, may be avoided by introducing scalar- matter couplings and/or 
by suitably choosing values of the slope parameters a and /3 that are phenomenologically 
viable. Extension to a more general situation with V{ip) V{(p, a) where the latter is a sum 
of exponential terms can be made rigourously by considering another scalar field a. This will 
appear elsewhere. 
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A. Appendix: 

In the appendices below we will present some technical details relevant to solving the field 
equations, along with a linear stability analysis of a critical solution. 

A.l Basic equations 

We find it convenient to define 

l^H^i^)^u{ip), ^ = ^'^x, (A.l) 

where a prime denotes the derivative with respect to the number of e-folds, N = j H dt = 
ln(a(t)/ao). A simple calculation shows that 

f = {uxy-e{ux), (p = {x' + ex)H^ ^y^ = ^-t^. (A.2) 
Equations (2.4)-( |2^ ) may be expressed as the following set of autonomous equations: 

= -3 + ^ + K^y + 2,K^xu + 317m + 317^, (A. 3) 

= 2e + 3 - K^{xu)' - K^{e + 2)XU + X^ - K^y + 2,Wm^rn + ^r, (A.4) 

= -fK^xx' + K^{£ + 3)72;^ + K^{y' + 2ey) + 3k^(1 + e)xu - S-qQ^mX, (A.5) 

where rj = 1, —2 and respectively for ordinary matter (wm = 0), stiff matter (wm = 1) and 
radiation {wr = 1/3). The continuity equation ( ^^ ) may be given by 

= 0;, + 2enm + 3(1 + W.m)^m + SrjQQrnX, (A.6) 
= 17; + 2(e + 2)17^. (A.7) 

In the case of minimal coupling, so Q = 0, we pi, oc e"^^^'^'"''^^. This last expression implies 
that the energy density due to radiation {wh = 1/3) decreases faster (by a factor of e^) than 
due to the ordinary matter (w^ = 0). 

A.2 A special solution 

As a canonical choice, consider the potential V{(p) and the coupling /{f) as the simplest 
exponential functions of the field: 

V{ip) = Vi e"^^/^« , f{^) = fi e"^/^» . (A.8) 

These choices may be motivated in the heterotic string theory as the first term of the pertur- 
bative string expansion [15]. One may solve the equations (|2.4D-(^^) by expressing V{ip) and 



/((/?) as in (|A.8D , with n = —m = 2, and dropping the effects of matter fields. The simplest 
(and also perhaps the modest) way of solving field equations is to make a particular ansatz 
for the scale factor, as well as for the scalar field, namely 

^^^^ijo) ' = {h>0,t>0), (A.9) 

/t -t\'^ t -t 
a = ao[— , ip = ipQ\n— (/i< 0, t < or < t < too)- (A.IO) 
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We can then reduce the field equations in the following form: 

= 7(1 - ih)vi + 2V,ti - " (A.13) 

H 

By solving any two of these equations, we find 

2 {hh-l)^^lK^ + Qh\h-\) 
Viti = 



ipQ=ipQ, (A. 14) 



for ^ 7^ —1, while 



,2_ 2 ^,2_4 6/o 

for /i = — 1. By fixing the values of Vq and ti, one can fix the value of lpq, in terms of the 
expansion parameter h, which itself needs to be determined by observation. The equation of 
state parameter can be written as 

w = -l + |. (A.16) 

Of course, w; < — 1 is possible for /i < 0. In this case the Hubble rate would increase 
with proper time t, since H = —h/{t^ — t). The observational results, coming from Hubble 
parameter measurements using Hubble space telescope and luminosity measurements of Type 
la supernovae, appear to indicate the value [39] 

-1.62 <w < -0.74 (A.17) 

at 95% confidence level. One may think of this value of tf; as two separate regimes: 

-1.62 < -u; < -1, -1< -u; < -0.74 (A.18) 

which, for the above critical solution, require that h < —1.0753 or 2.5641 < h < oo. That 
means, even if 7 > 0, as is the case for a canonical f, then the Gauss-Bonnet gravity could 
support both the effective phantom and the quintessence phase of the late universe. This 
was also one of the main points of the Gauss-Bonnet dark energy proposal made by Nojiri 
et al in [8], where the Gauss-Bonnet potential Vgb = f{'p)T^GB important, although it 
contributes to the effective action only subdominantly. 

Although the special solution above avoids the initial singularity at t = 0, it develops a 
big-rip type singularity at an asymptotic future i = too- Such solutions are typically only the 
fixed point, or critical solutions, which in general are not applicable once the matter degrees 
of freedom are taken into account, both for minimal and non-minimal couplings of (p with 
matter. Nevertheless, it may be interesting to check stability of such solutions under linear 
(homogeneous) perturbations, when V{ip) 7^ 0, which was not done in [8]. 
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A. 3 A linear stability analysis 

Our analysis in this subsection provides some new insights into the nature of instabihty of the 
critical solution under a linear (homogeneous) perturbation; the earlier analysis in [8] would 
be a special case of our treatment. 

With the ansatz ( A.9)-( A.IO), and the choice ( |A.8|) and m = —n = —2, the system of 
autonomous equations, in the absence of matter fields, is given by 

dx ^"^Lfox^ — 7x(2ux^ — ISuxifQ + Q(po) + 2(2?/ — 3ux) — 6n((/?o — 2xuipQ + x'^u) 
In ~ (po{2-f - 2-fxu + 3^2) ' 

du _ 2u [-^'^x'^(fo + u{2y - 3{ipo - x)u) + 27(2; - 2(pou)) 
m ~ ipo{2j - 2-fxu + 3u'^) ' 

where y = V/H"^. For the following special (critical) solution 

X = Xq = —, y = yQ = ——^ u = Uq = -i^ , (A.21) 

h 9304 

the right hand sides of the above equations are trivially satisfied, implying that x' = y' = 
z' = 0. 



(A.19) 
(A.20) 




Figure 20: Only for negative eigenvalues, satisfying the conditions A + D < and AD — BC > the 
linear perturbations become small and the system may be stable. 

Let us consider the following (homogeneous) perturbation around the critical solution 
( [Agp or ( [ATO| ): 

X = xq + 5x, u = uq + 5u. (A. 22) 

Since the functional form of the potential is already specified, the variation of y is encoded 
into the variation of 93. For stability of the solution, under a small perturbation about the 
critical or fixed point solution, requires 

A + D <0, AD-BC>0, (A.23) 
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where A, B, C, D are the eigenvalues of a 2 x 2 matric, as given by 



dN \dz J \Sz J ' \C D 

In the following, we define Vitl = k. We find 

^~ hE' ^~ ^E' ^~ 3hH3h^-3h-2k)E' ^~ hE ^^'^^^ 

where 

A = 18/i^(/i - l)(9/i2 - 4/i + 1) + 3/c(l - 24/i)/i^ + /i:2(132/i3 - 22h^ -7h + 1), 

(A.26) 

= 3{3h^ - 3/i - 2A;) (6/1^(3/1^ - 2h + 1) - k{22h^ - 13h + 1)) , (A.27) 
C = {36h^{3h^ - 1) - 2kh^{72h^ + 42/1^ - 12h{3h + 2) + 6) + 2k'^{2%h^ - 12h + 1)) 

x{9h^ -3h'^ -6hk + k), (A.28) 
D = 36h^{3h - l){h - if + 2k'^{24h^ - 68h'^ + 17h - 1) - 4:kh'^{36h^ - 93h'^ + 5lh - 6). 

(A.29) 

E = 18/1^(5/1^ - 4/1 + 1) + A;2(60/i2 - 16/i + 1) - Qkh^{2bh^ + 17/i - 2). (A. 30) 

We find the result in [8] by taking A: = In the case 7 > and < 0, all eigenvalues 
are negative, except when \h\ is small, \h\ < 2, or when the parameter k takes a large value 
(cf figure [20| ). While, for 7 < 0, two of them {B and C) take positive values, leading to a 
classical instability of the critical solution. Such a system is normally unstable also under 
inhomogeneous (cosmological) perturbations. In particular, for a large and positive potential, 
so that A; S> 0, only the /i > solution can be stable. 

A. 4 Fixed dilaton/modulus 

It is quite plausible that after inflation the field (p would remain almost frozen for a long time, 
so X ~ 0, leading to an universe filled dominantly with radiation and highly relativistic or 
stiff matter. Such a phase of the universe is described by the solution 

0^=1 - ^ o - -3 - 2pi e 3p2 e 

where Q,i = pi e^^ + p2 e^^ and pi are integration constants. One takes pi > 0, so that $7^ > 
and V{(j)) « A. The density fraction for a stiff fluid, f^^^^, decreases quickly as the universe 
further expands, ^> 0, thereby creating a vast amount of thermal radiation. 

While, at late-times, ip is again rolling slowly, such that x ~ 0. The cosmic expansion is 
then characterized by 

n-o-^ ^^3-45^ (A32) 



^The results here also correct errors/typos that appeared in the appendix of [8]. 
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where ^2 = e~ + S2 e and 5i are integration constants. Both 0,^ and fi^ decrease with 
the expansion of the universe, but unUke a universe dominated by highly relativistic fluids 
or stiff-matter, the radiation energy density decreases now faster as compared to the matter 
energy density, as it is red-shifted away by a factor of e~^. 
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